Slave spin cluster mean field theory away from half-filling: Application to the 
Hubbard and the extended Hubbard Model 
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A new slave-spin representation of fermion operators has recently been proposed for the half-filled 
Hubbard model. We show that with the addition of a gauge variable, the formalism can be extended 
to finite doping. The resulting spin problem can be solved using the cluster mean-field approxima- 
tion. This approximation takes short-range correlations into account by exact diagonalization on 
the cluster, whereas long-range correlations beyond the size of clusters are treated at the mean-field 
level. In the limit where the cluster has only one site and the interaction strength U is infinite, this 
approach reduces to the Gutzwiller approximation. There are some qualitative differences when the 
size of the cluster is finite. We first compute the critical U for the Mott transition as a function of 
a frustrating second-neighbor interaction on lattices relevant for various correlated systems, namely 
the cobaltites, the layered organic superconductors and the high-temperature superconductors. For 
the triangular lattice, we also study the extended Hubbard model with nearest-neighbor repulsion. 
In additionto a uniform metallic state, we find a y^(3) x y(3) charge density wave in a broad doping 
regime, including commensurate ones. We find that in the large U limit, intersite Coulomb repulsion 
V strongly suppresses the single-particle weight of the metallic state. 



I. INTRODUCTION 

The theoretical description of strongly correlated sys- 
tems, such as high temperature superconductivity, heavy 
fermions, and ultra cold atoms in optical lattices, etc., 
poses major challenges in field of the condensed mat- 
ter physics. These are all systems where the strength 
of the electron-electron interaction is comparable to 
or greater than the kinetic energy of the electrons, 
i.e., any theory based on a perturbative expansion 
around the non interacting limit is at least question- 
able. The non perturbative nature of the problems adds 
extreme difficulty to theoretical tools describing these 
systems. In recent years, several radically new and re- 
liable non perturbative approaches to the problem of 
strong correlations have been developed such as Dy- 
namical Mean-Feild Theory(DMFT)±, Dynamical Clus- 
ter approximation 2 -, Cluster-DMFT— , Variational Clus- 
ter Approximation (VCA)^, Two-Particle Self-Consistent 
Approach (TPSC)£; these new approaches have led to 
substancial progress in our understanding of these sys- 
tems. 

Some other non-perturbative semi-analytic approaches 
based on the idea of slave- variable representations of cor- 
related fermions have also been devised and have been 
used for decades now, in order to perform non-trivial 
approximations on many-body models. In this respect 
slave-bosons have been particularly successful. Their 
formulation in the limit of infinite correlation between 
the electrons^ can be systematically introduced as a sad- 
dle point approximation plus corrections, and has lead 
to much insight in the physics of the strongly corre- 
lated systems, most notably of heavy fermions. The 
alternative formulation that can treat finite interaction 
strength! cannot be controlled as a saddle point, but 
it turns out to be a very practical implementation of 



the Gutzwiller approximation. It has been generalized 
to many-orbital models^ and succeeded in capturing the 
essential of quasiparticle physics stemming out from the 
competition between interactions and derealization en- 
ergy. High energy features can be also studied from fluc- 
tuations around this mean-field. 

The main limitation of this last formulation is the fact 
that the number of slave-variables increases exponen- 
tially with the number of degrees of freedom in the mean 
field, making multi-orbital or cluster mean-field quickly 
intractable. 

A different approximation, based on quantum rotors as 
slave-variables^ has been devised that is much more eco- 
nomical since it introduces only one slave variable per 
site, dual to the total on-site charge. Still this tech- 
nique can only be used correctly at half-filling and can- 
not address orbital-dependent obscrvables or magnetic 
properties of the system. It has been nevertheless suc- 
cessfully applied to cluster men-fields recently^. Also an 
extension of this technique controlled by large degener- 
acy limits has revealed itself very powerful as an impurity 
solver—. 

Recently, a new representation of fermion operators 
that instead uses quantum spins as slave variables was 
proposed to study the multi-band Hubbard model at half- 
fillingi 2 -. In this paper, we generalize this representation 
away from half-filling and apply it to study Mott tran- 
sition on the different lattices and the charge denstity 
wave (CDW) transition on the triangular lattice. The 
Hubbard model plays the role of a standard model for 
correlated fermions on a lattice; it contains the band 
kinetic energy and the local on-site interaction. In or- 
der to study CDW, the Hubbard model was extended 
to include an intersite electron-electron interaction (V). 
This leads to the so called extended Hubbard model 
(EHM). Recently, EHM and its variant on the triangular 
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lattice have been extensively studied in the context of 
cobaltatealii^iS^iLia^a. 

The Hamiltonian for the extended Hubbard model 
(EHM) on a two dimensional lattice with sites labeled 
by i is 

<ij> i i 

+ V ^2 (m- 1)^-1), (l.i) 

<»i> 

where /x, i, [/ and V are the chemical potential, the 
nearest-neighbor hopping amplitude, the on-site local 
interaction U and nearest-neighbor interaction V, re- 
spectively, din {d\ a ) destroys (creates) an electron on 
site i with spin a, < i,j > denotes that the sum is 
over nearest neighbors only and the number operator is 

In order to treat this problem in the simplest approx- 
imation that is capable to lead to insight on the physics 
of short-range correlations, we employ a cluster mean- 
field approximation based on the slave-spin representa- 
tion. We have recently shown that cluster mean-field ap- 
proximation for bosons successfully describe the super- 
solid phase and phase diagram of bosons on triangular 
lattice^. 

In the following section, we introduce the method. In 
particular we introduce the gauge needed to its exten- 
sion off filled regimes. Sec. Ill presents the results on the 
Hubbard model, and Sec. IV those on the Extended Hub- 
bard model. We then summarize and conclude. Appen- 
dices contain various technical details such as the choice 
of gauge and the infinite U limit. 

II. SLAVE-SPIN MEAN FIELD THEORY 

Slave-spin mean-field theory^ is the ideal bridge be- 
tween the slave- variable techniques mentioned in the in- 
troduction, when taken at the mean-field level, in that it 
provides full insight in multi-obital and cluster cases, but 
still remains the most economical way to do this, since it 
introduces only one slave variable (a spin-1/2) for every 
degree of freedom in the mean-field cluster. In practice, 
for a single-site mean field of a one-band model, two slave- 
spins (one for spin-up electrons and one for spin-down 
electrons) are used, whereas for an N-orbital local mean- 
field or a N-site cluster mean-field of a one-band model 
the number raises only to 27V. Each slave spin increases 
the size of the Hilbert space by a factor of two. The 
gain is thus enormous compared with slave-boson rep- 
resentations because the number of bosons there grows 
exponetially. 

Where detailed comparison have been performed one 
finds, as discussed below, that the slave-spin mean-field 
reproduces the results of the Gutzwillcr approximation, 
even if a precise mapping has not yet been rigorously 
derived. 



A. Slave-spin representation for arbitrary filling 

In the slave-spin representation, we map the original 
local Hilbert space of the problem onto a larger local 
Hilbert space that contains as many fermionic degrees 
of freedom (named fi a ) as the original plus the same 
number of spin-1/2 quantum variables, one for each jiJ^. 
We then associate to every state of the original physical 
space one of the states in this larger space by using the 
correspondence: 

\ni = 1) \nl = 1, SI = +1/2), (2.1) 

\nf a = 0) ^ \nl = 0, SI = -1/2). (2.2) 

In words, when a local orbital and spin state is occu- 
pied then the corresponding slave-spin is "up" and if it is 
empty the slave-spin is " down" . With these one-particle 
states one construct the many-particle states as usual. 

The enlarged local Hilbert space contains also unphys- 
ical states such as \n[ a — 0,Sf a = +1/2) and \n{ a — 
1,S^ = —1/2). These unphysical states are excluded if 
the following local constraint is enforced at each site and 
for each a: 

flfi* = Sl+ 1 -. (2.3) 

We then have to map the operators onto operators that 
act in the enlarged Hilbert space. The electron number 
operator is easily represented by the auxilary fermions 
number, i.e., nf a = n{ a , but also by the z component of 
the slave-spin nf a = S 7 ( a + 1/2, thanks to the constraint. 
This allows us to rewrite the density-density interaction 
terms in the hamiltonian in terms of the spins only: 

Him = j EE 

i (7 

<i.j> o a 

For the non-diagonal operators we generalize the pre- 
scription of Ref.— , i.e. 

d la = f ia 2Sl, d\ a = fl2Sl (2.5) 

(where ] ia is the auxiliary fermion annihilation operator) 
to the more general one 

d itT = furOu,, 4 = flOj a (2.6) 

in which Oi a is a generic spin-1/2 operator, i.e. a 2 x 2 
complex matrix. 

Indeed it is easy to determine that the most general 
form for Oi a is 

Oia = ( J C f ) , (2.7) 
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where Ci a is an arbitrary complex number (When Ci a is 
not of unit modulus, there is no problem with anticom- 
mutation relations, if they are taken between physical 
states.), in order for the operator (|2.6p to have, in the 
physical states of enlarged Hilbert space, the same effect 
as the fermionic operators in the original Hilbert space, 
i.e.; 

d ta \ni = 0) = 0, d ia \ni = 1) = \ni = 0) 
4>t = 1) = 0, 4>t = 0) = \ni = 1) (2.8) 

The arbitrariness of the complex number Ci a is a gauge 
of our formulation and stems out from the fact that dif- 
ferent operators can have the same effect in the physical 
subspacc of the enlarged Hilbert space, while acting dif- 
ferently on the unphysical states. This difference does 
not have any effect as long as the constraint is treated 
exactly. In practice the local constraints are enforced via 
Lagrange multipliers and approximations have to be per- 
formed on these and on the Hamiltonian in order to solve 
the model. In these approximations the particular choice 
of gauge comes into play. Ci„ can indeed be tuned in order 
to give rise to the most physical approximation scheme, 
by imposing, for instance, that it correctly reproduces 
solvable limits of the problem, like the non-interacting 
limit. We will see that the correct choice of Ci a depends 
on the average occupation of the local state, and is such 
that it reduces to 1 at occupation 1/2, so that Oi a = 2Sf a 
and the prescription (|2.5[) used at half-filling in RefJ^ is 
correctly recovered. 

Finally, in the enlarged Hilbert space the Hamiltonian 
can be written exactly as: 

H = 1 E °l°j>rflfi° (2-9) 

<ij>a i<7 




E(E^) 2 + ^E(E^)(E^), 



i a <i,j> a a 

subject to the constraint (|2.3p . 



B. Mean- field approximation 

An approximation is now introduced, which consists in 
three main steps: 1) treating the constraint on average, 
using a static and site-dependent (but spin- independent, 
since we will not investigate here magnetic phases) La- 
grange multiplier A; 2) decoupling auxiliary fcrmions and 
slave-spin degrees of freedom and finally 3) treating the 
slave-spin Hamiltonian in a cluster mean-field approx- 
imation (CMFA), that takes into account the nearest 
neighbor correlations induced by V . 

After the first two steps, the total Hamiltonian can be 
written as the sum of the following two effective Hamil- 



tonians: 

H t =-t Qnflu* • //•'•• • - V; "< - 

<i,j>,a i 

(2.10) 

H s =- £ J ij OlO^+H.c. + Y J \{Sl + \) 

<r/>,cr i.a 

i <T <*,j> (7 a 

(2.11) 

The parameters Qij (effective hopping), Jy (slave-spin 
exchange constant) and Ai in these expression are deter- 
mined from the following coupled self-consistency equa- 
tions: 

Q l3 =< 0\ a O ]a > s , (2.12) 



Jij = * < flfja >/, (2.13) 



< nl >/=< SI > s +-, (2.14) 

where <>/, s indicates the effective Hamiltonian used for 
the calculation of the averages. We shall denote the near- 
est neighbor and next-nearest neighbor values of Qij as 
Q and Q' respectively. 

We are thus left with two coupled Hamiltonians: a 
renormalized free fcrmions Hamiltonian for the ji a and 
a lattice hamiltonian for the slave-spins that retains the 
full complexity of the original problem. We have thus to 
perform a further approximation, in this case the cluster 
mean-field on the spin Hamiltonian. 

A cluster with a finite number of sites only is consid- 
ered, within which interactions are treated exactly, and 
is embedded in the effective ("Weiss") field of its sur- 
roundings. A tiling of the original lattice is made, out of 
copies of the chosen cluster unit (Cluster shapes are cho- 
sen to respect lattice symmetry), assuming translational 
invariance in the supcrlatticc defined by this tiling, and 
this approximate Hamiltonian is used to calculate the 
mean-field average values. 

In practice this means that in this approximation an 
effective Hamiltonian of a finite cluster is enough to rep- 
resent the physics of the full lattice and that the "Weiss 
fields" are calculated using this same Hamiltonian (i.e. 
sclf-consistcntly) that represents also the surroundings 
of the cluster unit and not only the cluster unit itself. 

Mathematically, we consider the following Hamiltonian 
for the slave-spin cluster C : 

h° = H > [*. ft + E h ^°l + H - c - + E w > 

<ij>eC ieC \cr ie C 

(2.15) 
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FIG. 1: Modulus and phase of the gauge c that determines 
the choice of the proper hopping operators in the enlarged 
Hilbert space, in order for the CMFA to reproduce the non- 
interacting limit. 



where hi„ and hf are the effective fiefds of the sur- 
roundings, that are determined by the following self- 
consistency conditions: 



jn.n.i 



(2.16) 



(2.17) 



jn.n.% 



where the prime over the sum means that sites j inside the 
cluster are excluded. We solve the spin Hamiltonian on 
cluster size N c = 3 for the triangular lattice and N c = 4 
for the square lattice. 

It is useful to underline the role of two key quantities, 
in characterizing the physics of the system. It can be 
shown that Z =< O la > 2 is the quasiparticlc weight, 
while the effective mass enhancement is set by the ef- 
fective hopping renormalization Qij. The two quantities 
coincide if the mean-field approximation on the slave spin 
hamiltonian is taken at the single-site level. Thus, they 
both vanish in the Mott insulating phase. This amounts 
to neglecting all number fluctuations within the Mott 
phase. This is too crude of an approximation especially 
when close to the Mott transition. On the contrary in 
the CMFA that we consider here these two quantities are 
distinct and one can have e.g. a Mott transition where 
the mass stays finite as we will see in the following. 

We use Z as an order parameter: Z ^ indicates 
a metallic state, while Mott/CDW insulating behaviour 
corresponds to Z = 0. 



FIG. 2: Illustration of the lattices with hopping amplitude t 
and t' . (a) Vo x \/3 sublattice decomposition of the anistropic 
triangular lattice (ATL) (b) and (c) two sublattice decompo- 
sition of the anistropic frustrated square lattice (AFSL) and 
the isotropic frustrated square lattice (IFSL). A, B, and C 
indicate the sublattice decomposition. 



C. Choice of the gauge Ci a 

We now discuss how to fix the gauge represented by 
the complex number ci a . 

The physical condition that we choose to impose is 
that our CMFA reproduces correctly the non-interacting 
limit, i.e. when U = V = 0, 



Qi 



z = i. 



(2.18) 



for any given filling , so that Ci a = c(n ia ). 

In the single-site approximation c can be chosen purely 
real and it can be determined analytically (as detailed in 
Appendix A). It takes the form: 



1 



y/n{l - n) 



1 



(2.19) 



More generally c has to be determined numerically by 
solving the mean-field equations at U = V = and im- 
posing the conditions (|2.18p and is a complex number, 
i.e. c = |c|e i<6 In Fig. Q]we show \c\ and as a function 
of for a triangular cluster and for the single-site result, 
both on a triangular lattice. 

We note that in both cases at half-filling a = 1 and 
= 0, and Oia coincides with the form chosen in Ref^, 
as anticipated. 



III. HUBBARD MODEL 

The Mott transition, i.e, the metal-insulator tranisition 
driven by the strength of electron-electron interaction in 
a homogenous phase, has been studied in a great detail 
using various approaches such as slave bosons, DMFT 
and its extensions. In this section, we revisit the Mott 
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transition on the lattices shown in Fig.l. The control 
parameters are interaction strength U/t and frustration 
strength t' /t, the ratio of next nearest neighbor to near- 
est neighbor hopping amplitude. As a function of these 
parameter, the Hubbard Model at half filling has, within 
CMFA, four possible phases : a paramagnetic metallic 
phase, a paramagnetic insulating phase, and insulating 
antifcrromagnctic phase, and (in the presence of frustra- 
tion ) an itinerant antifcrromagnctic phase. However, we 
shall be concerned here with the transition between the 
paramagnetic metal to parmagnctic insulator. We study 
the paramagnetic solution by enforcing the spin symme- 
try hence avoiding the opening of a full spectral gap due 
entirely to magnetic ordering. 

The single-site cluster N c =l mean- field theory of 
slave spin representation gives the same results of the 
Gutzwiller approximation. In this regards CMFA pro- 
vides a way to go beyond the Gutzwiller approximation. 

First, we discuss the Mott transition on the isotropic 
triangular lattice. To get the uniform phase solution, 
we enforce the Lagrange multiplier Ai and complex num- 
ber c to be the same for every site within the cluster. 
In FigO we plot Z and Q as a function of U at x=0 
for cluster sizes N c = 1,3. The critical value U/t, at 
which the Mott insulating phase occurs is 16.2, in the 
single-site (JV C = 1) approximation, while it is around 
15.1 in the three sites (JV C = 3) CMFA. The short range 
correlations, which are built in the CMFA, supress the 
critical value U by 6%. The critical value of U obtained 
from other methods such as DMFT-exact diagonalisa- 
tion (8 site)22, exact diagonalisation calculation for 12 
site clusters^, and cluster-DMFT(CDMFT)^ are 15, 12 
and 10.5 respectively. In CDMFT the tranisition is first 
order. For N c = 1, the slave spin approach is identical 
to the Gutzwiller approximation^, Q and Z are identical 
and they should vanish at the same critical value of U. 
We show, for N c = 3, Q as a function of U. It can be 
seen that it continues to be non-zero in the Mott insu- 
lating phase, and behaves as a t/U, as expected from the 
fact that the average kinetic energy is non-zero in a Mott 
insulator. 

It should be noted there is a substancial difference 
between U c obtained from 3-site slave-spin CMFA and 
3- site CDMFT. It is because of CDMFT captures the 
fermionic quantum dynamics more accurately. 

We now examine Z as a function of dopping in the limit 
U — » oo since this quantity can be obtained in closed form 
in the Gutzwiller approximation. We show in Appendix 
B that for cluster size A c =l, one recovers precisely the 
Gutzwiller approximation result Z = 2x/(l + x). In 
CMFA, we can ask how Z is modified in the presence 
of short-range correlation effects. It can be seen in Fig0] 
that the short range correlation effect on Z is apprecia- 
ble for moderate to large doping x and enhances Z in 
comparison to cluster size N c = 1. 

We now move on to the the dependence of the Mott 
transition on lattice and frustration. In the absence of 
magnetic frustration on a bipartite lattice, one expects 




FIG. 3: The order parameter Z and the effective hopping Q 
as function of U. U is measured in unit of t. 
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FIG. 4: The order paramter Z as a function of x in the large 
U limit for N c =l,3. 



to find an antiferromagnetic ground state at low temper- 
ature. Ideally, the Mott transition can occur in system 
where antifcrromagnctic correlations are frustrated. In 
the t — t' Hubbard model on the square lattice, a next- 
nearest neighbor hopping t' frustrates antifcrromagnctic 
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FIG. 5: The order parameter Z as function of U at t! — 0.4 
for the lattices ATL, AFSL, and IFSL. U is measured in units 
oft. 



5 10 , 15 20 25 



0.2 0.4 , 0.6 0.8 1 



FIG. 6: The order parameter Z as a function of U for various 
values of t' for IFSL. U is measured in units of t. 



FIG. 8: The order parameter Z and the effective hopping i 
as function of U. U and t' are measured in units of t. 




FIG. 7: The nearest neighbor effective hopping Q as function 
of U at t' = 0.4 for AFSL and IFSL lattices. U is measured 
in units of t. Inset shows the next nearest neighbor Q' as a 
function U in unit of if at t = 2.5. 



correlations. By studying the lattices shown in Fig. 2, 
we thus investigate the effects of frustration on the Mott 
transition in the half-filled t — t' Hubbard model. For 
t' = the lattices shown in Fig.l correspond to the un- 
frustratcd systems and the effect of the frustation can be 
systematically studied as t' is increased to its maximal 
value t' = t. FigjS] displays the order parameter Z or 
single particle weight as a function of U at t' = for 
various lattices. At this value of t' ', the critical value U 
for the Mott transition is the lowest on the anistropic 
triangular lattice (ATL), while it is the highest on the 
isotropic frustrated square lattice (IFSL). As the frustra- 
tion t' increases, the critical value of the Mott transition 
increases as shown in Figl6] for the isotropic frustrated 
square lattice. This increase with t'/t is also seen in 
the Variational Cluster Approximation^. In Fig[7l we 
show the nearest and the next-nearest neighbor effective 
hopping Q and Q' of auxilary fermions. It can be seen 
that deep in the insulating phase they behave ast/U and 
t' /U respectively. Non-zero values of Q and Q' in the 
insulating phase signal that auxilary electrons (not the 
physical electrons) have a Fermi surface (with Luttinger 
Volume). It also implies, in contrast with infinite dimcn- 
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FIG. 9: This phase diagram displays critical values of U and 
V where Z vanishes at half-filling. U and V are measured in 
unit of t. 



sion (where single-site mean-field theory is exact), that 
in finite dimension the effective mass does not diverge in 
the insulating phase, despite the fact that Z — > 0. 

Finally, we show in Fig. 8 the phase diagram in U — t' 
plane for the above mentioned three lattices. One notices 
that the maximally frustated lattice, the triangular lat- 
tice, has the lowest critical value of the Mott transition, 
while isotropic frustrated square lattice has the highest 
critical U. 



IV. THE EXTENDED HUBBARD MODEL 

A. Uniform Phase 

In this section, we consider the extended Hubbard 
model on the isotropic triangular lattice for different 
ranges of parameters U, V , and doping (x) in the uni- 
form phase, by enforcing that the Lagrange multiplier be 
the same at every site on the cluster and thus avoiding 
charge ordering. Let us first cxamincthc combined effect 
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FIG. 10: Z as a function of dopping x at U = 100£ for var- 
ious values of V. W is the full bandwidth of the isotropic 
triangular lattice (W=9t). 

of U and V on Z at x = 0. For given U or V, wc compute 
the crititcal value of U or V at which the Mott transition 
occurs. This study leads to the uniform ground state 
phase diagram in the U — V plane that is shown in Fig. 
[5J For U < 7.5 the system is in the metallic state for any 
values of V. For 7.5 < U < 15, the system enters into the 
Mott insulating phase upon increasing V. We note how- 
ever that there is a 'reentrant' structure of the metallic 
phase at still larger V. This 'reentrant' structure emerges 
when U and V are comparable. It is because of V com- 
pensates the effect of U. And moving a nearest-neighbor 
to have a doubly occupied site, as in a metallic phase, 
may become energetically favourable since the repulsion 
on the nearest-neighbor is comparable to that on-site. 

From study of Sodium Cobalt Oxide in Refi2S, it ap- 
pears that there is a large supression of the valence- 
band width— by an order of magnitude compared with 
the local density approximation (LDA) band structure 
calculation^!. Ref.— suggested that such large renormal- 
ization of the hopping may be caused by V, and thus Q 
was studied as a function of x for different values of V 
by means of the Jastrow-Gutzwiller (JG) wave function. 
The Q in JG wave function study is equivalent to Z in 
our case. Slave-spin CMFA should be more accurate than 
the Jastrow-Gutzwiller approximation since, it captures 
the short-range correlation effect of V in a better way 
because this term is treated exactly on the cluster. 

In FigfTUl we show Z as a function of x for different 
values of V at (U = 100). It should be noted that the 
value of Z vanishes at the commensurate dopings x — 1/3 
and x = 2/3 when V takes its largest value, V = 7W/4 
(where W = 9t is the full bandwidth of the isotropic 
triangular lattice). At doping 2/3 the dominant configu- 
rations at large V on any triangle are (|, j, f) (J., f, |) 
(I, I, I). Now Z involves Slipping a spin. So we have to 
make transition to states like ((J., !, J.) or (f, f, |) etc. 
These have a higher energy in the presence of V. Similar 
argucments holds for at doping 1/3. Z vanishes in our 
case around V = 7W/4, which is quite a large value in 
comparsion to the JG study, where it occurs at V ~ W. 



This implies that JG study overestimtes the effect of the 
short-range correlation of V. 

The effect of V on the effective hopping tQ is not as 
strong as we observe on Z (not shown) since on a sin- 
gle triangle, there is no cost to move the particle via a 
kinetic move, eg: Q, j, f) — > (J., |, j). It is true on a 
single triangle not connected on anything else, but the 
mean-fields connected to the triangle will have small ef- 
fect, which manifests itself by a small supression in the 
effective hopping. Going beyond a single triangle how- 
ever, it is clear that the effective hopping (bandwidth) 
will be suppressed, a possiblity which wc do not consider 
here. 



B. The CDW instability 

In this section, we study the instability of a non- 
ordered phase toward a CDW in the presence of V. We 
determine for a few dopings the ground state phase di- 
agram in the U — V plane of the system that has the 
V3x V3 ordering pattern. We solve the above mentioned 
equations allowing for site-dependent Lagrange multipli- 
ers. The procedure is as follows: We allow the complex 
number gauge c and the Lagrange multiplier A different 
for each sublattice and Q's and J's are different for every 
bond on the cluster. We fit the filling dependence of the 
magnitude of the complex number c and its phase to ob- 
tain c for different sublattices, based on the c(n) relation 
obtained on the noninteracting system.. 

Fig.(|lip shows the resulting CDW phase diagram in 
the U—V plane for three special values of x. The transi- 
tion from metallic to CDW phase is first order. We also 
note that the effective mass 1/Q diverges at the transi- 
tion. We note that the lowest value of V c is at x = 1/3. 
We also do not find that dopings x = 1/3 and 2/3 are 
playing any special role, as was suggested in the JG study. 
It should also be noted that we find the CDW state at 
x = 0.5(not shown) in contrast with the prediction of 
the uniform phase in our study (where Z never vanishes) 
and in the JG studyi^. The slave-boson mean field study 
of Re£i£ also predicts a phase diagram similar to ours. 
However, our method captures the short-range effect of 
correlations in a better way than Re£i£. We suspect that 
our phase diagram does not match with the diagram pro- 
posed by the JG study because it underestimates the ef- 
fect of V. 



V. CONCLUSION 

We have presented an extension of the slave-spin for- 
malism away from half-filling by introducing a gauge vari- 
able. And we have shown how to solve the resulting 
model in the cluster mean-field approximation (CMFA). 
While in the single-site mean-field approximation the 
gauge variable can be chosen as pure real number, it is a 
complex number in the cluster approximation. This num- 
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FIG. 11: Phase diagram in U- V plane for x = 1/3, 1/2, 2/3. 
The CDW phase is above the lines for the corresponding fill- 
ings. 



be generalized to tackle the t-J Model, and other strongly 
correlated models. 
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APPENDIX A: CHOICE OF THE GAUGE C IN 
THE SINGLE-SITE MEAN-FIELD 



ber changes from a pure real to a pure imaginary number 
as we move from the haff-filled to the empty lattice. The 
advantage of this method lies in the fact that the short- 
range correlations can be properly taken into account. In 
the single-site approximation for the Hubbard model, we 
found analytically that in the infinite U limit, the single- 
particle weight Z reproduces the Gutzwiller result. In 
the CMFA, short-range correlations modify this result. 
The modifications are more important for intermediate 
dopings but they are never very large. 

We have applied this approach to the Hubbard and to 
the extended Hubbard Model. In the case of the half- 
filled Hubbard Model, we have revisited the Mott transi- 
tion on three class of lattices: anistropic triangular lattice 
(ATL), High-T c lattice(ISFL), and organic superconduc- 
tor lattice (AFSL). We have done a detailed study of the 
critical value U where the Mott transition occurs as a 
function of the frustration strength t' , and have shown 
that the effect of t' in the presence of the short range 
correlations is to increase the critical value for the Mott 
transition U c {t'). 

We have also studied the extended Hubbard model in 
two dimensions in the uniform phase and shown that 
there is a reentrant structure between the insulating and 
metallic phases when U and V are comparable. We have 
also shown that dopings 1 /3 and 2 /3 play a special role in 
the uniform phase. The quasiparticle weight can vanish 
at these dopings. 

For the extended Hubbard model, we have found two 
ground state phases on the triangular lattice: the metallic 
and the V3 x V3 CDW state in a broad doping regime. 
At the present level of approximation, we found that, 
contrary to the uniform phase, dopings 1/3 and 2/3 in 
the CDW state do not play a special role. 

Finally, we point out that this method can be used 
to study magnetic phases. That has been left for future 
work. It can also be applied to study the physics of the 
multiband Hubbard model away from half-filling and can 



In the single-site approximation, we can determine the 
gauge c analytically^. 

The non-interacting single-site slave spin Hamiltonian 
H s reads^ 



H s = hO^ + h*0 + X(S Z + -), 



(Al) 



where O is defined as in cq. (|2.7|) . The single-site fermion 
part of the Hamiltonian is simply spinlcss non- interacting 
fermions. The physical spin index a is suprcsscd in H s 
since for U = upspin and downspin fermions are de- 
coupled, so that we can diagonalize the hamiltonian for 
one slave-spin in the S z — ±1/2 basis. The ground state 
eigenvalue ecs and the corresponding eigenstate are 
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z 2 l = -R 



\GS) 



JV 
— a* 

N 



(A2) 



(A3) 



with N = y2R(q + R) and a = h + ch*. 

The expectation value of S z and O in the ground state 



(S, 



_\_ 

4R 



and 



(O) 



ca 



2R 



(A4) 



(A5) 



The Lagrange multiplier depends on the density n and is 
adjusted in order to satisfy the constraint equation: 



n~2 = {Sz 



A 

AR 



(A6) 
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We want to tune c in order to match the condition that 
in the limit U — the rcnormalization factor Z must be 
unity: 



where j indicates the neighbor of site i. Diagonalizing 
H s , we obtain the ground state eigenvalue and eigenvec- 
tor, i.e. 



Z =< O > 2 = 



4R 2 



= 1 



(A7) 



We can easily eliminate A from these two conditions, 
by squaring eq. (|A6[) . We are left with the following 
expression for c: 

- n 2 (A8) 



lea* + a\ 2 



If we choose c to be real then h and a are also real. 
Then, the expression for c in the closed form is 



1 



y/n(l-n) 



- 1. 



(A9) 



Note that this result is independent of h. 

This cannot be done in the cluster case, since also the 
condition Q=l has to be imposed and c has to be chosen 
complex in order to satisfy this further equation. 



with 



A 1 



\GS) 



N 




G.S 



2\a\ 



Hence we can determine (S*) and (O'). 



(Ot) 



2 N 2 
a*e-G.s 
e 2 GS + 2\a\ 2 



(1 + c). 



(B4) 



(B5) 



(B6) 



(B7) 



APPENDIX B: SLAVE-SPIN FORMULATION OF 
THE INFINITE- [/ LIMIT OF THE HUBBARD 
MODEL 

We derive here the analytic expression for Z as a func- 
tion of doping in the infinite U limit and in the single-site 
approximation. 

In this limit, no double occupancy is allowed so that 
the interaction term is replaced by a projector that en- 
forces this constraint. In order to do this we replace Oi a , 
as defined as in eq. (|2.7|) . by 



O ia = {\-Sl s )O l 



(Bl) 



where a = — a. We thus obtain, for the single-site mean- 
field spin Hamiltonian H s : 

H s = J2h*di+H.c. + \J2(S* + b (B2) 



with 



■J2 J ij(Oj.) (B3) 



The Lagrange multiplier is fixed by the constraint equa- 
tion that depends on the chosen filling: 



1 



2 e 2 GS + 2\a\ 2 



(B8) 



We can calculate the renormalization factor Z 
using 



|(Ot)|2 



Z = \(0)\ 2 = \l + c\ 2e f^n 2 
\a\ 2 

= |1 + c\ 2 n{l - 2n) 



(B9) 



Using the one-band prescription (see appendix [XJ c = 
. 1 — 1, we then obtain 



1 - 2n 2x 



1 — n 1 + X 



(BIO) 



(where x is the total doping, 2n = 1 — x). That is pre- 
cisely the result of the Gutzwiller approximation. 
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